Phonon decoherence of quantum entanglement: Robust and 

fragile states 



T. Yu and J. H. Eberly 

Rochester Theory Center for Optical Science and Engineering 
and the Department of Physics & Astronomy 
University of Rochester, 
Rochester, NY 14627, USA 

(Dated: February 1, 2008) 

Abstract 

We study the robustness and fragility of entanglement of open quantum systems in some exactly 
solvable models in which the decoherence is caused by a pure dephasing process. In particular, for 
the toy models presented in this paper, we identify two different time scales, one is responsible for 
local dephasing, while the other is for entanglement decay. For a class of fragile entangled states 
defined in this paper, we find that the entanglement of two qubits, as measured by concurrence, 
decays faster asymptotically than the quantum dephasing of an individual qubit. 

PACS numbers: 03.65.Yz, 03. 67. -a 
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I. INTRODUCTION 



Central to quantum information and quantum computation is the concept of entangle- 
ment of qubits [1] . In ideal situations, entangled quantum states would not decohere during 
processing and transmission of quantum information. However, real quantum systems will 
inevitably be influenced by surrounding environments. The interaction between the envi- 
ronment and a qubit system of interest can lead to decoherence. This is manifest in the 
loss of unitary evolution [2, 3]. The decoherence process varies for different quantum states. 
On the one hand, the pointer basis is formed by those states that are unaffected by the 
environmental variables, so in this sense they constitute a set of robust states [2]. If such 
robust states are also entangled, the entanglement is expected to be stable [4, 5, 6, 7, 8]. 

Our intuition strongly suggests that a specified entanglement, as a nonlocal property of 
a composed quantum system, should be very fragile under the influence of the environment. 
This fragility is a main obstacle for the realization of practical quantum computers. Among 
the various proposals to combat the decoherence in quantum computing and quantum infor- 
mation processing are ion traps and nuclear magnetic resonance (NMR) [9, 10, 11, 12, 13]. 

The main purpose of the present paper is to focus on the issue of fragility. We show 
in a specific case that the environment affects the type of coherence called entanglement 
quantitatively more severely than it affects the coherence associated with off-diagonal matrix 
elements of a single qubit. We show that this demonstration can be carried out without 
approximation for a non-trivial open system described by a reasonably general Hamiltonian. 
Here we write i^tot ^ the sum of Hamiltonians for the system itself, the environment, and 
the coupling between them (we use h — 1): 

Htot ^H + J2 {9\La\ + glL'^ay) + ^ u;a4«a, (1) 

A A 

where H is the Hamiltonian of the system of interest. The coupling operator L is a system 
operator coupled to the environment, and the gx are coupling constants. 

We cannot discuss all issues in full generality. Specifically, we present a toy model in 
which the environment is represented as a pure dephasing process. Our final conclusions, 
based on calculations of concurrence [14] , will apply to a two-qubit pair, and it is fair to say 
that entanglement of two qubits is a well-discussed topic [15]. However, we think that the 
way in which entanglement itself decays (or doesn't decay) when a system is exposed to a 



2 



nonlocal noisy channel is worth a quantitative examination. 

The organization of the paper is as follows. In Sec. II, we present an exactly solvable 
model and the solution of the exact non-Markovian master equation. In Sec. Ill and IV, 

we first present a two-qubit model and then provide some detailed studies of robust and 
fragile entangled states that are initially pure, as they relax to mixed states. We compare 
the entanglement decay time with the local dephasing time in Sec. V. We conclude in Sec 
VI. 

II. AN EXACTLY SOLVABLE MODEL 

Equation (1) describes a standard model for open quantum systems in the system- plus- 
rescrvoir framework in which a system is coupled linearly to a bath of harmonic oscillators, 
the excitations of which we can interpret as phonons or photons, for example. In any event, 
the bath has distributed eigenfrequencies ujx and creation and annihilation operators a\, ax 
satisfying [0^,0]^/] = ^a,a'- We assume that the system and the environment are initially 
uncorrelated: p(0) = Ps(0) ® Pbathi^), where Pbathi^) represents the thermal state of the heat 
bath at temperature T. 

In the following, we will consider a specialized model such that H = W and L = L'^ 
and the two self-adjoint operators satisfy [L, H] = inl, where k is a constant, and / is the 
identity operator acting on the Hilbert space of the system. A particular interesting case is 
when K — O. That is, H and L commute with each other. Let us note that, except for these 
conditions, for the time being, we do not assign operators H and L any concrete forms. 

The exact non-Markovian master equation for the model presented above takes the fol- 
lowing form [16]: 

= -i[H, pt] + F{t)[Lpt, L] + F*{t)[L, ptL] (2) 
+G(t)[pt,L] + G*{t)[L,pt]. 

where 

F(t) ^ a(t,s)ds, and G(t) ^ k [\(t, s)(t - s)ds, (3) 
Jo Jo 

where a{t, s) — r]{t, s) -\- ^{t, s) is the bath correlation function at temperature T: 

V{t,s) = l^|^Arcoth(^^^^cos[a;A(i-s)], (4) 
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u{t,s) = -J2\gx\'^sm[u;xit- s)]. (5) 

A 

Note first that this master equation is in a time-local form. So finite memory and thus 
non-Markovian effects are encoded in the time-dependent coefficients F{t) and G{t) [17, 18]. 
In particular, the G{t) term is a pure non-Markovian term which does not appear in the 
Markov approximation, when a{t — s) — >■ 6{t — s). This term gives a non-Markovian phase 
shift. 

Let us now turn to the special case: k — 0, i.e., the system's Hamiltonian H and the 
coupling operator L commute with each other. Suppose the states \n) are simultaneous 
eigenkets of H and L: 

H\n) — En\n), and L\n) = /„|n). (6) 

We can explicitly solve the corresponding non-Markovian master equation (2), and the 
solution is surprisingly simple: 



Pnm{t) = {n\pt\m) (7) 

where the functions Fii{t) and Fj{t) are the real and imaginary parts of F{t), respectively. 
The first exponential factor in (7) represents a phase shift and the second one introduces 
decay, i.e., the decoherence effect. In what follows we always assume that F{t) has an 
asymptotically positive real part ensuring that the decoherence is an irreversible process. 

We can see from the solution (7) that eigenvectors of L are robust states. Precisely, for 
any initial pure state of the system \ip{0)) = |n), we have p(0) = \tjj{0)){^jj{0)\ = p{t) = 

m))m)\- 

III. TWO-QUBIT SYSTEM 

So far, we have not made any concrete assumptions about the structure of Hamiltonian 
H and the coupling operator L. In order to discuss entangled states, we have to specify the 
Hamiltonian. This section is devoted to discuss a simple yet interesting example. 
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The system we consider consists of two coupled qubits A and B, where the Hamiltonian 
for two qubits is taken to be nonhnear and nonlocal: 

H = ujAcrf + ujbct^ + J<J 
and the coupling operator is given by 



(8) 



(9) 



The coupling operator L commutes with the Hamiltonian H and this guarantees that the 
energy is conserved at any time. Thus the decoherence is a pure dephasing process in which 
the loss of quantum phase of the system into the environment is the source of decoherence. 
The dephasing is a key issue in practical implementations of quantum computers [11, 12]. 

As shown in the last section, the eigenvectors \n) of L are robust states. The solution of 
(2) for the two-qubit system is given by 



Pt 



Pii(O) Piz{t) Pu{t) 

P2l{t) P22(0) /923(0) P24{t) 

Psiit) P32(0) p33(0) p34{t) 

Pii{t) Pi2{t) Pisii) Pu{0) 
where we have employed the "standard" eigenbasis: 

|1)ab = |++), |2)as=| + -), 

|3)aB-|-+), H)aB=\--): 

and where the matrix elements are: 

xe-^/o^«(^)<^Vi2(0), 

^-i{Ei-E3)t-i4. J*F,{s)ds 
^-i{E2-E4)t+i4 J*Fi{s)ds 

xe-^/o^«W^>24(0), 

g-i(E3--E4)t+j4 J^Fi{s)ds 

xe-^i'o^«(^)<^V34(0). 



Pisit) = 

Pu{t) = 

P24{t) = 

Pu{t) = 



(10) 



(11) 

(12) 
(13) 

(14) 

(15) 



Here the eigenvalues of H and L are given by 

Ei = UJi^- UJ2 + J, E2 = UJl - UJ2 - J, 

Es = —cui + 102 — J, £"4 = —CiJi — u;2 + J, 



and 



respectively. 



h — 2, h — h — 0, I4 — —2, 



IV. DEGREE OF ENTANGLEMENT: ROBUST VS FRAGILE 

In order to quantify the degree of entanglement, we will adopt the concurrence C defined 
by Wootters [14] . The concurrence varies from C = for an unentangled state to C = 1 for 
a maximally entangled state. For qubits, the concurrence may be calculated explicitly from 
the density matrix p for qubits A and B: 

C(p) = max{0, Ai - A2 - A3 - A4}, (16) 

where the quantities Aj are the square roots of the eigenvalues in decreasing order of the 
matrix 

^ = P«®<7>*(<7X), (17) 

where p* denotes the complex conjugation of p in the standard basis. 

The most general pure states in the case of the two-qubit model can be written as 



\^)ab = ai|l)As + a2\2)AB + a^\2,)AB + a4|4)AB, (18) 
where Yh=i — 1- The concurrence of the pure state (18) is simply given by [14] 

C(|*)ab) = 210203 -aia4|. (19) 

Thus, the pure state (18) is entangled if and only if 

0104 7^ 0203. (20) 

By a robust entangled state we mean one whose entanglement will not decay to zero in 
temporal evolution (10). In what follows, we will prove that the following two special cases 
with 04 = and Oi = 0: 
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\iPi)ab = ai\l)AB + a2\2)AB + a3\3)AB, 
\i^2)AB = 02|2)aB + a3|3)AB + a4|4)AB, 



(21) 
(22) 



are robust entangled states if 0203 7^ 0. In doing so, let us compute the concurrence of the 
density matrix with the initial state (21). The density matrix for qubits A and B at time t 
is given by 

I flip Pl2{t) pisit) 
P2l{t) |«2P a2«3 

psi(t) aa*^/? I flap 


It is easy to check that the concurrence 



Pt 



(23) 



C{pt) = C(po) = 2|a2a3| 



(24) 



That is, state \iI'i)ab has robust entanglement (which is more than simply being a robust 
state). Similarly, we can show, for the entangled pure initial state (22), that the degree of 
entanglement is completely preserved. 

The above result is slightly surprising since the pure states (21) and (22) are the su- 
perpositions of the eigenvectors of L with different eigenvalues. One would expect that 
decoherence would eventually degrade the degree of entanglement. In fact, if we use 



P{t) = Tip', 



(25) 



to quantify the loss of purity of a quantum state, for the initial pure state (18) we have: 



4 r ft 

P{t)^ E laiHa/exp -{h - Ijf dsFnis) 

i,j=l 



(26) 



As a measure of purity of a quantum state, it is easy to see that < P < 1 where P = 1 if 
and only if p represents a pure state. Prom (24) and (26) we can easily see that in temporal 
evolution the purity of the state (21) deteriorates, but the degree of entanglement remains 

constant. 

In what follows, we will investigate another class of entangled pure states whose en- 
tanglement tends to vanish under the influence of the environment. Specifically, we will 
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demonstrate, for the entangled bipartite pure states (18) with — and 02 = 0: 

= ai\l)AB + a2\2)AB + a4\4:)AB, (27) 

\(f>2)AB = ai|l)AB + 03|3)AB + a4|4)AB, (28) 

that the entanglement will vanish after an entanglement decay time denoted by r^. We refer 
those states as fragile entangled states. 

Now we are in the position to discuss the entanglement decay of the fragile states by 
explicitly computing the concurrence. The density matrix with the initial entangled state 
(27) at t is given by 

|aiP Piiit) puit) 
p2i{t) |a2p P24(t) 



Pt = 



P4l{t) p42{t) 1 04 1 2 

and the concurrence can easily be obtained 



(29) 



C{pt) = 2\pu{t)\ = 2|aia4|e-^^/o 



(30) 



The time scale for the entanglement of a fragile state decaying to zero is determined by the 
function Fji{t). For large times, the function Fji{t) T = r]{t)dt. Then in this long 
time limit, the entanglement decay time can be identified as 



r-' = ler. 



(31) 



V. ENTANGLEMENT DECOHERENCE VS LOCAL DEPHASING 

The dcphasing rate of an individual qubit can directly be estimated from the density 
matrix for qubit A, which can be obtained from the density matrix (10) by further tracing 
out the variables of qubit B, and vice versa; that is, p'^ — Tr^p, p^ — Tr^p. The reduced 
density matrix for qubit A is thus obtained from (10): 

Pll + P22 Pn + P24 
P31 + P42 P33 + P44 

Thus, the dephasing rate denoted by for qubit A is determined by the off-diagonal elements 
of pt, 

|pi'2l-|Pi3 + P24|~e-^/o'^^^«(^). (33) 



pt= 
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Similar analysis applies to qubit B as well. Clearly, the dephasing time depends on the 
behavior of the function FR{t). Similar to the entanglement decay time, by ignoring the 
details of the heat bath we can immediately identify the dephasing time t<^ for the large 
times as 



We have thus used the exphcit solution of master equation (10) to evaluate the time 
development of both the degree of entanglement and the dephasing rate. As seen from 
the expressions (30) and (33), the entanglement for qubits A and B and local quantum 
coherence for an individual qubit, say qubit A, decay at different rates. We have shown, for 
a large class of fragile entangled states, that the entanglement decay time is shorter than 
the local dephasing time on which the quantum coherence of each local qubit is destroyed. 
Moreover, it can easily be shown, in the all cases with the entangled initial states (18), that 
the entanglement decay time is not longer than the dephasing time, i.e., 



VI. CONCLUSION 

Robust and fragile states have been discussed in some exactly solvable models. Par- 
ticularly, for the two-qubit toy model, we have investigated the dynamics of the robust 
and fragile entangled states in terms of a measure of entanglement called concurrence. We 
have emphasized and identified that there exist two different decoherence time scales — 
entanglement decay time and local dephasing time. For the coupling considered here, we 
find that the entanglement decay occurs faster than the local dephasing for a large class 
of fragile entangled states. If the size of the active qubits were to greatly increase, the en- 
tanglement decoherence time would be expected to become exceedingly small, reflecting the 
classical limit [19]. Our example supports one's intuition about the non-local nature of the 
entanglement. 

Since the amount of entanglement contained in an entangled quantum state is dependent 
on the choice of a specific measure of entanglement, the entanglement decoherence rate 
is also dependent on such a choice. We also emphasize that classification of both robust 
and fragile entangled states will depend on the concrete form of the interaction between 



r^' = 4r. 



(34) 



(35) 
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qubits and the environment. In addition, we do not expect that our toy model exhibits 
all interesting aspects of the decay processes of the quantum entanglement. We do beheve, 
however, that the fast decay rate of quantum entanglement is a generic feature in a variety 
of physical processes where decoherence is important. 

Finally, it may be worth noting that our two-qubit model allows two competing pro- 
cesses: creation and annihilation of entanglement [20, 21]. Thus, the maximal entanglement 
generation under the dephasing process is an interesting problem that will be addressed in 
future publications. 
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